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Abstract. This paper deals with the existence of traveling fronts for the 
reaction-diffusion equation: 

^-Au = h(u,y) t e R, x = {x u y) e R N . 

We first consider the case h(u,y) = f(u) — ag(y)u where / is of KPP or bistable 
type and limij / i_>_|_ 00 g(y) = +00. This equation comes from a model in population 
dynamics in which there is spatial spreading as well as genetic mutation of a quan- 
titative genetic trait that has a locally preferred value. The goal is to understand 
spreading and invasions in this heterogeneous context. We prove the existence of 
threshold value ao and of a nonzero asymptotic profile (a stationary limiting so- 
lution) V{y) if and only if a < ao- When this condition is met, we prove the 
existence of a traveling front. This allows us to completely identify the behavior of 
the solution of the parabolic problem in the KPP case. 

We also study here the case where h(y, u) = f(u) for \y\ < L\ and h(y, u) m —au 
for \y\ > L2 > L\. This equation provides a general framework for a model of 
cortical spreading depressions in the brain. We prove the existence of traveling 
front if L\ is large enough and the non-existence if L2 is too small. 



1 Introduction 

This paper deals with the existence of bounded traveling fronts for the reaction- 
diffusion equation 

^-Au = h{y,u) t eM, x= (x u y) eR N . (1) 



1 



The function h will be of three different forms in this paper. The first two concern 
non-linear terms h(y, u) = f(u) — ag(y)u where / : M — > R is C 1 , and either of KPP 

type, or of bistable type and g : R" -1 — > M + is C°, 5(0) = and g ^^ + °°> +00. 
The existence of traveling front depends on the value of a > 0. The third case we 
consider here is when h(y,u) = f(u) for \y\ < L\ and h(y,u) < —mu for \y\ > L2 
where 0<Li<^2<oo are given parameters and / is of bistable form and 
h(y, u) + mu — > for \y\ — > +00. We study the existence of traveling fronts depend- 
ing on the value of L\ and Li. 

The problems we study in this paper bear some similarities with the question 
of traveling fronts in cylinders of [3J. However there are important differences that 
have to do with the fact that the cross section in [3j was bounded and only the Neu- 
mann condition was considered there. Whereas here, the problem is posed in the 
whole space and the solution vanish at infinity in directions orthogonal to the direc- 
tion of propagation. We follow the same general scheme as in [3J and in particular 
the sliding method. But some new ideas are also required. In particular, first, we 
treat directly the KPP case without the approximation of the KPP non-linearity by 
a combustion non-linearity as in [3j . Then in the approach of Bcrestycki - Nirenbcrg 
[3] to traveling fronts in cylinders for the bistable case, a useful result of H. Matano 
|27j was involved in the proof. Here, we rely on stability ideas but also use energy 
minimization properties to bound the speed of the solution in the finite domain 
approximation. In particular, we do not use the precise exponential behavior that 
was used in [3j. Actually the developments of this method that we present in this 
paper can be used to somewhat simplify parts of |3j. They can also be applied to 
traveling fronts in cylinder with Robin or Dirichlet boundary conditions, lj 

Equation (JT|) in the first case comes from a model in population dynamics [20) 
that we briefly describe now. Let u(t, x, v) represent the density of individuals at 
time t and position x that possess some given quantitative genetic trait represented 
by a continuous variable v G K. The latter could be for example the size of wings or 
the height of an individual. We assume that individuals follow a brownian motion 
(i.e. they diffuse) in space with a constant diffusion coefficient v, reproduce identi- 
cally and disappear with a growth rate k(x, v) that depends on the position x and 
on the trait v. Furthermore, they also reproduce with mutation that is represented 
by a kernel K(x,v,w) and disappear due to competition with a constant L > 0. 
Thus, one is led to the following equation for u: 

dtu(t, x, v) — vA x u(t, x, v) = k(x, v) u(t, x, v) , . 

+ J K(x,v,w) u(t, x,w) dw — u(t, x,v) J Lu(t, x,w) dw. 

We assume moreover that there exists a most adapted trait </> = <f>(x) that may 
depends on the location x. The farther the trait of an individual is from the most 
adapted trait, the larger the probability of dying and not reproducing. Thus the 
growth rate can be written for example as k(x, v) = a — b \v — (j>(x)\ 2 with a and 
b > 0. Non-local reaction-diffusion equations of this type are quite difficult to handle 
from a mathematical standpoint as shown in |I2j where behaviors very different from 

lr The construction of traveling fronts for Neumann and Dirichlet conditions in cylinders given 
by [41| appears to be incomplete. Indeed, the continuity of the function <f> on page 515 is not 
established before using Dini's Theorem to derive Lemma 3.2 there. 
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those in local equation are brought to light. In a forthcoming numerical study [5], we 
show that depending on the value of a, the traveling fronts may not be monotonous 
anymore. This fact was also established in a different context for some related types 
of non-local equations in |12j . This leads us to introduce a simplified version of this 
model that emphasizes propagation guided by the environment. First, we assume 
that mutations are due to a diffusion process represented by a Brownian motion in 
the space of trait v. Furthermore, we assume that </> is linear. Then a rotation in 
the variables (xi,y) allows one to reduce the problem to the case where the most 
adapted trait is y = 0. Therefore we assume <j>(x) = and © can be rewritten as 

dtu(t, x, v) — vA x . v u(t, x, v) = (a — b\v \ 2 )u(t, x, v) — u(t, x, v) / Lu(t,x,w) dw. (3) 

J w 

Lastly we assume that competition is only between individuals sharing the same 
trait which leads us to equation 

dtu — vA XlV u = (a — Lu)u — 6|u| 2 w. (4) 

Equation (TTJ is a generalization of this equation. In [50], the authors observe nu- 
merically a generalized transition front spreading along the graph of <f> for equation 
([2) (see [SI M EH HH] or [37] for the definition of generalized transition fronts). Here 
we want to prove theoretically (i) that there exists such a front for equation (fTJ at 
least for some values of the parameter a > and (ii) that extinction occurs if a 
is too large. The latter condition can be interpreted as saying that the "area" of 
adapted traits is too thin compared to the diffusion. To remain consistent with the 
biological motivation, we only consider here non-negative and bounded solutions of 
©■ 

Other types of models related to this one have been proposed in the literature. 
For example, the model developed by Kirkpatrick and Barton in 1997 [25] also 
studies the evolution of a population and of its mean trait. The main difference is 
that they have a system in u and v where u represents the population and v the 
mean trait is described by a specific equation. This model has been further explored 
many times since [22] [24] . It is worth noting that these models use the same type of 
non-linearity for the adaptation to the environment and model the mutation with 
the Laplace operator as well rather than the integral operators. 

This type of reaction-diffusion process in heterogeneous media also arises in 
many contexts in medicine. An important class of such models was treated in 
|15[ 133] . They deal with the propagation of a cortical spreading depression (CSD) 
in the human brain. These CSD's are transient depolarizations of the brain that 
slowly propagate in the cortex of several animal species after a stroke, a head in- 
jury, or seizures [3S]. They also are suspected of being responsible for the aura in 
migraines with aura. CSD's are the subject of intensive research in biology since 
experiments blocking them during strokes in rodents have produced very promising 
results [121 130] • These observations however have not been confirmed in humans and 
the existence of CSD's in the human brain is still a matter of debate [2"§1 12"31 [Tl I3T)] . 
Since very few experiments and measurements on human brain are available be 
it for ethical or for technical reasons, mathematical models of a CSD is help in 
understanding their existence and conditions for propagation. In such a problem, 
the morphology of the brain and thus the geometry of the domain where CSD's 
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propagate, is believed to play an important role. 



The brain is composed of gray matter where neuron's soma are and of the white 
matter where only axons are to be found. The rodent brain (on which most of the 
biological experiments are done) is rather smooth and composed almost entirely of 
gray matter. On the opposite, the human brain is very tortuous. The gray matter 
is a thin layer at the periphery of the brain with much thickness variations and 
convolutions, the rest of the brain being composed of white matter. According to 
mathematical models of CSDs [TTJ [3H1 [Ml SO] , the depolarization amplitude follows 
a reaction-diffusion process of bistable type in the gray matter of the brain while 
it diffuses and is absorbed in the white matter of the brain. The modeling of CSD 
hence leads one to the study of equations of the following type: 



Here, / is of bistable type and \y\ = L corresponds to the transition from gray mat- 
ter to white matter. In [15], this equation was studied to prove that the thinness of 
the human gray matter (L small) may prevent the creation or the propagation of 
CSDs on large distances. It was proved by studying the energy in a traveling refer- 
ential of the solution of §5$ with a specific initial condition. The special case of (O 
for N — 2 was described more completely in [TH]. In [33]; a numerical study shows 
that the convolutions of the brain have also a strong influence on the propagation 
of CSD. Finally, in [16], the effect of rapid variations of thickness of the gray matter 
was studied. 

Finally, let us note that the same kind of equation arises in the modeling of 
tumor cords but with a slightly more complicate KPP non-linearity. We plan to 
investigate this model in our forthcoming work [6]. 

As already mentioned, the study of propagation of fronts and spreading proper- 
ties in heterogeneous media is of intense current interest. For instance, the existence 
of fronts propagating in non- homogeneous geometries with obstacles has been estab- 
lished in Berestycki, Hamel and Matano [7]. Definitions of generalized waves have 
been given by Berestycki and Hamel in [9] and |10j where they are called general- 
ized transition waves. Somewhat different approaches to generalizing the notions 
of traveling fronts have been proposed by H. Matano [28] and W. Shen [37] . The 
existence of fronts for non- homogeneous equations are established in [3T] and [12"] . 

Let us first introduce some notations before stating the main results. 

Notations. We note x = (xx,y) € where x% € M and y £ R^" 1 . Hence x 
is the space variable in M. N , x\ is its first coordinate and y is the vector of IR^ -1 
composed of all the other coordinates of x. As usual Br = B(0, R) denotes a ball 
of radius R centered at 0, but here it will always mean the ball in R N_1 . 

First we are interested in solutions of 



— - Au = f(u)l\ y \ <L - aul ]y] > L 



t G R, x = (x x ,y) € R- 



(5) 




— - Au = f(u) - ag{y)u, x = {x\,y) <E Mr, (eR 
u > 0, u bounded, 



(G) 
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with a > 0. We will assume that / : R — > R is C 1 and satisfies either one of the 
following conditions: 

/(0) = /(l) = 0, / > on (0, 1) and f(s) < f'(0)s for a G (0, 1), 

or 

there exists 9 e]0, 1[ such that /(0) = /(0) = /(l) = 0, 
/ < on (0, 6) and / > on (0, 1) and /'(0) > 0, f (1) > 0, 

Jo > °- 

The first case will be referred to as the KPP case and the second one will be called 
bistable case. Since we are only interested in solutions of (O in [0, 1], we will further 
assume that f(s) < for s > 1. Moreover we always assume 

g : R^ 1 -> K+ is continuous, g(0) = 0, g > on R w_1 \ {0} (7) 

and 

lim g{y) = +oo. (8) 

|y |— >+oo 

Taking g(y) = \y\ 2 and f(s) = os(l — s) yields the particular case of equation (QJ. 

This paper is concerned with the long term behavior of ([6]) and with the existence 
of curved traveling fronts, i.e. solutions u(t, x) = U{x\ —ct, y) with c£la constant 
and U : M. N — > R such that the limits lim s _ i .± 00 U(s, .) exist uniformly and are not 
equal. Regarding these fronts, our main results are the following. 

Theorem 1.1. If f is of KPP type, there exists ckq > such that: 

• For a > ckq, there exists no traveling front solution of , 

• For a < ao there exists a threshold c* > such that there exists a traveling 
front of speed c of equation ^ if and only if c > c* . 

This existence theorem gives us information on the behavior of the solution of 
the parabolic problem. In this paper we prove the following theorem: 

Theorem 1.2. If f is of KPP type, for uq £ L°° , there exists a unique solution 
u(t,x) of 

d t u — Au = f(u) — ag(y)u on (0, +oo) x R N , 
u(0,x) = uo(x) on M. N . 

• If ot > ao, it verifies u(t,x) t ~ >+ °°^ q uniformly with respect to x £ R w . 

• If a < ao and uq € C^R^) with uq < V where V = V(y) is the unique 
positive asymptotic profile (stationary solution ). 

for any c > c* lim sup u(t, x) = 0, 

|xi|>ct 

for any c with < c < c* lim sup \u(t, x) — V(y)\ = 0. 

t ^+ 00 |x 1 |<ct 



This means that there is a threshold value «o such that for a > ao, there is 
extinction. On the contrary, when a < ao, there is spreading and the state V(y) in- 
vades the whole space. The asymptotic speed of spreading is then c* . The property 
of asymptotic spreading is in the same spirit of the theorem of asymptotic speed of 
spreading in cylinders established by Mallordy and Roqucjoffre in |26j . 

Regarding the case of bistable / we have the following result: 
Theorem 1.3. /// is of bistable type, there exist a* > a* > such that 

• For a > a* , there exists no traveling front solution of (J5]) , 

• For a < a*, under condition \40\ of Section^ there exists a traveling front u 
of speed c > solution of ((6]) . 

Lastly, the model for CSD's leads one to equations of the type 

d t u- Au = h(y,u) x = (xi,y)eR N . (9) 

where h{y, u) that verifies 

h(y,u) = f(u) for \y\ < Li 
h(y,u) < —mu for \y\ > 

|y| — ^-)-00 1 

h{y, u) + mu > uniformly for u £ R 

where < L\ < Li < 00 and m > are given parameters and / is of bistable form. 

In this paper we prove the following Theorem. 

Theorem 1.4. There exist critical radii < < L* < 00 with the following 
properties: 

• For L12 < L*, there is no traveling front solution of ©• 

• For L\ > L* (independently of L2 ), assuming that there is a unique stable 
asymptotic profile of (f50|) . there exists a traveling front of speed c > solution 
of ©. 

This result completes the study in [T3] on the existence of CSD in the human 
brain. Indeed in |15j the transition from gray to white matter was instantaneous 
when biologically there is a smooth transition from gray to white matter. This The- 
orem confirms the intuition that CSD's can be found in part of the human brain 
where the gray matter is sufficiently thick but they can not propagate over large 
distances due to a thin gray matter in many parts of the human brain. 

The paper is organized as follows. In section[5]we state some preliminary results 
that will be used in the sequel. Section [3] is dedicated to the study of the existence 
and uniqueness of non-zero asymptotic profiles for a traveling front solution of © . 
In section [S]we study the large time behavior. There we prove extinction if a > a a 
and convergence towards the front of minimal speed if a < a . Then, section [S] is 
devoted to the study of the asymptotic profiles in the bistable case and section [7] to 
the existence of traveling front for a < a* in the bistable case. Lastly, in section [8] 
wc describe the precise problem arising in the modeling of CSD's and state our 
main result in this framework. 
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2 Preliminary results 

In our proofs, we will several times need the exponential decay of the asymptotic 
profile which can be easily proved from the following theorem established in [14] . 

Theorem 2.1. Let v G iJ ; 2 oc (IR^ ) be a positive function. Assume that there exists 
7 > and C > such that 

Vz G R N , v(x) < Ce^M and liminf > 7- 

|:c|->oo v(x) 

Then, lim v(x)e^ VA = 0. 

¥00 

This result is established in [14], lemma 2.2. In the context of equation ((T|), we 
thus have the following corollary. 

Corollary 1. Let u be a non-negative and bounded solution of 

Av + f(v) - agv = on R N_1 . 

Then, for any 7 > there exists C > such that 

< v(y) < Ce- 7|y| and \Vv(y)\ < Ce'^ 1 . 

Proof. The estimate on v comes directly from Theorem 12.11 and the estimate on 
|Vii| derives from standard global LP estimates. □ 

3 The case of a KPP non-linearity. Asymptotic 
profiles. 

In this section, we are interested in the asymptotic profiles of a traveling front 
solution of © as x\ — > ±00. Hence, we are looking for solutions of the following 
equation 

AV + f(V)-ag{y)V = 0, y eR N -\ 
V>0, Abounded. 

We assume that / : E ->• R is C\ 

/(0) = /(l) = 0, / > on (0, 1) (11) 

and 

f(s) < f'(0)s for s G (0, 1) and s € (0, 1] ^ LU. decreasing. (12) 

s 

Since the constant function is always solution, the problem is to know when 
there exist non-zero solutions. As we will see here, the existence of such a positive 
asymptotic profile is characterized by the sign of the principal eigenvalue of the 
linearized operator around 0. We now make precise this notion. 
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3.1 Principal eigenvalue of the linearized operator 

To start with, let us define the natural weighted space 

n = {ve H 1 ^- 1 ) , e £ 2 (R Ar_1 )} 

and its associated norm. For v e %, we set \\v\\u = {\\v\\ 2 H1 + Hv^llia)^- The 
linearized operator about is Lip = —Aip + (ag{y) — f'(0)*)<p for <p G H. We are 
interested in the eigenvalues of L. Even though the problem is set on all of l* -1 , 
the term in ag(y) yields compactness of the injection % L 2 (R N ^ 1 ). Hence the 
existence of a principal eigenvalue is obtained as usual. 

Theorem 3.1. Let us define 



The operator L has a smallest eigenvalue 

X a = inf R a (<p). (13) 

<p£H\{0} 

Moreover there exists a unique positive eigenfunction associated to X a of L 2 -norm 
equal to 1, called tp a in the following,. The eigenspace associated to X a is spanned 
by ip a . 

Proof. The proof is classical due to the compactness of TL ■— >• L 2 (R N " 1 ). We refer 
for example to (3T]. □ 

Remark 1. If g(y) = \y\ 2 , the problem can be rescaled and we obtain the harmonic 
oscillator for which principal eigenvalue and eigenfunction are well known [35j . In 

that case, we have X a = (N - 1)^/5 - /'(0) and ip a = (^f) e~^l 2 . 

Since the existence of a positive solution of (|10[) will depend on the sign of the 
principal eigenvalue, the following proposition describes the behavior of X a as a 
function of a. 

Proposition 1. The function a > X a is continuous, increasing and concave for 
a e (0, +oo). Moreover lim Q _>o X a = —f'(0) and for a large enough X a > 0. 

Proof. Let us fix a > and 77 > 0. Equation (|13l) shows that 

A Q+)) < J |V^ Q | 2 + {(a + rf)g - /'(0))^ 2 =X a +r,J g(y) V 2 a . 

Similarly, we obtain X a < X a+V — ij J gip 2 a+n - From this we derive that 

< rj J g<p 2 a+1) < X a+V - X a < ?] J gip 2 a . 

This and similar computation for X a — A a _ l; yields that a i— > X a is increasing and 
locally Lipschitz on (0, +00). 
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Concavity is classical. It suffices to observe that for each fixed ip, 

is an affine function of a and that A Q = inf^g^^o} R a ((p). 

In order to prove that A Q > — /'(0), for any e > choose a function ^> e of 
compact support with HV'elU 2 = 1 an d J \S7tp E \ 2 < e. Let supple C Br f . From 
JT3J we get 

-/'(0) < A Q <e + amax.g-/'(0) 

So for any a < 



max BH( g 

-/'(0)<Aa<-/'(0)+2e. 

Now we claim that Aq, > for large enough a. Argue by contradiction and 
assume that A Q < for all a € (0, +oo). Since 

0>A Q = J |V^ a | 2 +a fgt&-f'(0), 

therefore 

g^i < -no) 

a 

and (p a — > in L 2 (M JV ~ 1 \ Br) for all i? > 0. Furthermore, ip a is bounded in "H and 
up to extraction we can assume that ip a converges strongly in L 2 (R W_1 ). thus <p a 
converges to in L 2 but it is impossible since J (p a 2 = 1 for any a > 0. □ 

Corollary 2. There exists ao > such that A Q < for a < cxq, X ao = and 
A a > for a > ao- 

3.2 If g vanishes on B ro 

The main part of the proof is still correct if g vanishes on B ro but the result is 
slightly modified. 

In this section, we assume that there exists ro > such that ([7]) is substituted 
by the following assumption 

g : l^ 1 -> R + e C°, g = on B ro and g > on R^ 1 \ B ro . (14) 

We define Aa the principal eigenvalue of the Laplacian on B ro with Dirichlct bound- 
ary conditions, i.e. 

J -A</>o = A A 0o on B ro , 
1 0o = on dB ro . 

In this case, the principal eigenvalue of the linearized operator about is well 
defined and Proposition Q] becomes 

Proposition 2. The function a t— > A Q is continuous, increasing and concave for 
a g (0, +oo), and lim Q _>o A Q = — /'(0). -/Vow there are two cases: 
i) If /'(0) < Aa, £/ien /or a /arge enough A Q > 0. 
#7'(0) > A A , then X a < /or aZZ a > 0. 
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Proof. The proof of the first part of the proposition is exactly the same as in Propo- 
sition [T] We just have to prove i) and ii). 

i) We assume that /'(0) < Aa and argue by contradiction that A Q < for any 
a € (0, +00). As in the proof of proposition [TJ we have 

/ < i/'(0) 

but this yields ip a — > in L 2 (Jst N ~ 1 \ Br) for a — > +00 and any R> r . 

As before ip a is bounded in % and up to extraction, we have A Q — > A < and 
the weak convergence in T-L and strong convergence in L 2 of ip a toward 4>. The limit 
4> verifies / </> 2 = 1, 4> = for \y\ > ro and 

-A<f>-f'(0)<l> = \<t> 

Thus (f> G Hl(B ro ) verifies -A</> = (/'(0) + A)</> on B ro and since A < 0, it is 
impossible. 

ii) By taking ip = 4>q in the Rayleigh quotient p^|) . where </>o is the principal 
eigenvalue of the above problem in B ro with Dirichlet boundary conditions, we see 
that A Q < A A - /'(0) < for all a > 0. □ 

In the following, we will not state the results specifically for this case (fT4)) . 
However, the proofs and results developed here carry over to this case with the 
obvious modifications. 



3.3 Existence of non-zero asymptotic profile 

Theorem 3.2. For a > ao, there is no solution of (|10|) . where ao is defined in 
corollary \^ For a < ao, there exists a unique positive solution of (|10l) . 

Proof. Let us fix a > ao- Then A Q > 0. Assume by contradiction that there exists 
a solution V of ([10]) . Then the strong maximum principle shows that V > 0. 

Since tp a is an eigenfunction of the linearized operator L and V is solution of 
([TO]), we have 

J (AV + f(V) - agV) <p a = 

= J (A<p a + (/'(0) - ag)ip a + \ a ip a )V 

Now from corollary [U V and W are rapidly decreasing for \y\ — > 00 and so we 
can apply Stokes formula J AV (f a = fVAip a . It yields J (f(V) — f(Q)V)ip a — 

A Q J tp a V but f(V) - f'(0)V < since / is of KPP type and \ a > thus a 
contradiction is obtained. 



We now turn to the case a < ao- For a < ao, the eigenvalue A Q is negative. 
Setting V_ = sip a with e > 0, we get 

AV + f(V) - agV = -X a eip a + f(e<p a ) - f\G)e<p a > 
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if e > is chosen small enough. Hence V_ is a sub-solution of ([TU]) . The constant 
function 1 is a super-solution and V_ < 1 if e is small enough. Therefore by the sub- 
and super-solution method, there exists a solution V such that < V_ < V < 1. 

Now consider V and two non-zero solutions of (fTUfl . We argue by contradic- 
tion and assume that V ^ W. Then for example fi = {y € K^ -1 , U(y) < W(y)} is 
not empty. Introduce a cutoff function /3 S C°°(R) with /3 = on (-co, 1/2], /3 = 1 
on [1, +oo) and < < 4 on (1/2, 1) and for all e > 0, let us set /3 s (s) = [3 (f ). 
Using equation (fT0|) . we have 

J (—VAW + AVW)f3 s (W - V) = J (Vf(W)- f(V)W)p e (W -V) 



> / (Vf(W) - f(V)W) 



£2 



by Lebesgue's dominated convergence theorem. Owing to corollary [T] V, VV, 
and VVF have exponential decay and thus Stokes formula can be applied and we 
obtain 

J (—VAW + AVW)/3 E (W — V)=J fi' E (W - V)V(W - V). (WW - WW) 

f3' e (W - V)V\V(W - V)\ 2 - I p' e (w - V)(W - V)V(W - v).vv . 



In the term li the integrand satisfies 

\P' E (W - V)(W - V)V(W - V).VW\ < 4|V(W - U)|.|VW| 

Therefore by Lebesgue's Theorem of dominated convergence, we infer that I2 — > 0. 
Next the term 1\ satisfies 1\ > 0. Consequently, we may write: 

> J (vf(w) - wf(v)) = jT UQp- Hp-^j VW 

which is a contradiction in view of (|12p as W > V in 57. Hence V — W and the 
non-zero solution is unique. □ 

The last point concerns the stability of the asymptotic profiles for a < cto. Let 
us start by studying the energy of V. For hi £ M, we define the energy 

J a (w)= f \Vw(y)\ 2 + ^g(y)w 2 (y)-F(w(y))dy (15) 
where F(u) = £ f(t)dt. 

Theorem 3.3. For a < ao, the unique positive solution of (1101) V is stable in the 
energy sense, i.e. V is the global minimum of J a and, furthermore J a (V) < = 
J a (0). 

Proof. Owing to the maximum principle, solutions of (fl"U|) are between and 1. 
Hence we can modify / on ] — 00, 0[ such that it becomes odd and as a consequence, 
F can be considered as even. Since A Q the principal eigenvalue of the linearized 
operator about the zero solution is negative for a < ao, cannot be the global 
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minimum of J a . Now J a admits a global minimum that will be called V for the 
argument. One can prove that \V\ is also a global minimum of J a and hence \V\ is a 
positive solution of (|10j) . By uniqueness, \V\ = V and thus V is a global minimum 
of J a . Since is not a global minimum, necessarily J a (V) < = J a (0). □ 

We now conclude with the linearized stability of V. 

Theorem 3.4. For a < ao, consider the linearized operator about V and denote 
Ai[V] the principal eigenvalue of this operator. Then Ai[V] > 0. 

Proof. Denote by tp a positive eigenfunction associated with Ai [V] and assume by 
contradiction that Ai[V] < 0. If M[V] < 0, it is easy to see that for e > small 
enough V + etp < 1 is a sub-solution of (fTU|) . From there, it would follow that there 
exists a solution of (fT0|) between V + etp and 1 but this contradicts the uniqueness 
of V. 

Now if Xi[V] — 0. With the same notations as above, tp and <j> verify 

-Aip = f(V)ip - ag(y)ip (16) 

and 

- Acp = (/'(0) + \ a )<p - ag(y)cp (17) 
Multiplying ([T6]) by <p and (fT7|) by tp and integrating by parts yields 

J f'(v)n = J (f(0) + \ a )tp<p. 

Since <ptp > 0, we obtain f(V(y)) = /' (0) + A Q for all y g l"" 1 and this is 
impossible since A a < and V{y) —> for \y\ — > +oo. 

□ 



4 Traveling fronts for a KPP non-linearity 

This section is devoted to the definition of a speed c* for which a traveling front 
of equation © exists for a £ (0, ao)- The threshold of existence of the non-zero 
asymptotic profile is called ao as in the previous section. For < a < ao, V 
denotes the unique non-zero asymptotic profile. As shown in the previous section, 
the energy of the non-zero profile J a (V) is negative. 

A curved traveling front of speed c is a function u(x\ —ct, y) solution of equation 
(|6|) and connecting the non-zero asymptotic state V to 0. Thus we are looking for 
a solution of 

!— Au — cd\u + ag(y)u = f(u), x = (xi, y) g 
u(x u .) Xl ^- QO ) v 7 u(x u .) Zl ~ >+0 °> uniformly in y g R^ 1 , (18) 
u > 0, u bounded 

where c g R is also an unknown of the problem. 

The construction of c* in Theorem 11.11 uses the sliding method following ideas 
of [3]. Note however that there are important differences with [3J. In that paper, 
the KPP case is derived by first solving the "combustion non-linearity" and then 
approach the KPP non-linearity as a limiting case of truncated functions. Contrary 
to [3J here, we derive directly the existence of a solution of the KPP case. Actually 
the method we present here can be applied to somewhat simplify the proof of [3J in 
the KPP case for cylinder with Neumann conditions. 
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4.1 Problem on a domain bounded in x±. 

Let us fix a > 1 and c € R for this subsection and consider the following problem: 

-Au — cd\u + ag(y)u = f(u), x = (xi, y) G (—a, a) x R-^ -1 
u(-a,-) = V, u(a,-)=0, (19) 
u > 0, it bounded. 

The aim of this subsection is to prove the following theorem: 

Theorem 4.1. There exists a unique solution of 119]) . denoted u c a in the following. 
This solution decreases in the X\-direction, i.e. d\u c a < 0. Thus < u c a < V 
on (—a, a) x R^ -1 . Moreover c u c a is decreasing and continuous from K to 
L°°{[-a,a] x R^" 1 ). 

To prove this theorem, we require the following two propositions. 

Proposition 3. Let u be a solution of (fT9]) . Then u(x\,y) < V(y) for {x\,y) £ 
[-a, a] x R^- 1 . 

Proof. Let M > 1 be such that u < M and consider tpR defined on Br the largest 
solution of 

f -Ayip R + ag(y)ip R = f(ip R ) for y e B R , 

\ip R = M iory€dB R , < ip R < M. l_J ' 

Here we think of / as having been extended by outside [0, 1]. Since f(s) < for 
all s > 1, we observe that: 

• by the strong maximum principle, < ip R < M on B R . 

• since V < 1 < M and V is a sub-solution of ([2U)) , through monotone iterations 
we have V < ip R . 

• if R' > R, ip R > is once again a sub-solution of on B R and thus ip R > < ip R 
on B R . 

• therefore ip R tends to a function when R — > +oo and through local elliptic 
estimates, this function is a non-zero solution (> V) of the asymptotic problem 

(flO|) . By uniqueness, we obtain tp R R ^ +co > y 
Now we consider the problem 

j-Aio - ediw + ag(y)w = f(w) for x <G (—a, a) x B R , 

1 w = M for x e (—a, a) x w = for xi = ±a, j/ G 



The solution u of (jT9J) is a sub-solution of ([21) and the constant function M is a 
super-solution. Using monotone iterations starting from the super-solution M, we 
build the same sequence as previously (for problem <\20]i ) since by induction the 
solutions do not depend on x\ G (—a, a). Hence the sequence converges toward ip R 
and we have u < ip R < M. Now letting R —> +oo yields u < V. □ 
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Proposition 4. Let R be such that g(y) > ^ for y ^ Br where K is the Lipschitz 
norm of f on [0,1]. We set f2 = / X (K^" 1 \ Br) where I is an open bounded 
interval of IR. 

Suppose u and v £ C 2 (f2) n C°(f2) are solutions of 

— Aw — cdiw + ag(y)w = f(w) on 57 (22) 

and u < v on dfl. Then u < v on O. 

Proof. By contradiction, suppose this is not true. Due to corollary[T]and proposition 
[3l u{x\,y) and v(x\,y) converge uniformly to for \y \ — > +oo. Consequently, there 
exist (xo, yo) S ^ such that 

> min(i) — u) = (v — u)(xo,yo). 
n 

Since (xo,yo) G ^ we have d±(v — u)(xo,yo) = and A(v — u)(xo,yo) > 0, and 
subtracting the equation ()22[) with u from the one with v, we obtain 

ug{\yo\){v - u)(x ,y ) > f(v(x ,y )) - f(u(x ,y )) > -K\(v - u)(x ,y )\ 

which is impossible since ag(|yo|) > K and (v — u)(xo,yo) < 0. □ 

Let us now turn to the proof of Theorem 14. II using sliding method. 

First u(x, y) = V(y) is a super-solution, is a sub-solution and < u, so by 
monotone iterations, there exists a solution u of (|19|) . 

Lemma 4.2. Assume u and v are two solutions of h!9\). Then 

v{x\ + h,y) < u(xi, y) for all h € [0, 2a) and all (x\, y) G [—a, a — h] x K.^ -1 . 

Proof of the lemma. By proposition [3l we have < u < V (resp. < v < V) and 
using the strong maximum principle, we obtain < u < V (resp. < v < V) on 
(-a, a) x R N -\ 

For h € [0, 2a), let I h = (—a, a — h) and for (xx, y) 6/|,x IR^ -1 , set Vh{x\,y) = 
v{xi + h, y). 

Let us fix R > such that g(y) > — for y £ Br. By compactness and continuity 
of u and v, there exists e > such that Vh < u on Ih x Br for any h such that 
2a — e < h < 2a. Proposition |4] shows that Vh < u on Ih x R w_1 for any h > 2a — e. 
This enables us to define 

fr* = inf{/i > 0, u fe < n on 4 x M^" 1 }. 

Let us prove that h* = and argue by contradiction that ft.* > 0. By continuity, 
v h , < u onlftTx l w_1 . 

Suppose that min u — Vh* > 0. This would imply that for h* — h > small, 

I h , xB R 

min u — Vh > and by Proposition [H v/j < u on //j x in contradiction 

Ih x Br 

with the definition of h* . Therefore min u — Vh* =0. This implies the existence 

I h , xB R 

of (x\,y*) € Ih- x such that Vh* [x\ , y* ) = u(a;|,y*) (note that w — Vh* > 
for x\ = —a or a — h*). Writing in the usual way that u — Vh* is solution of a 
linear elliptic equation in Ih* x R^ -1 and u — Vh* > with m — Vh* vanishing at 
the point (x1,y*), the strong maximum principle implies that u — Vh* = which is 
impossible. □ 
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Applying the preceding lemma with h = yields the uniqueness of the solution 
of (fP9]l . Taking u = v = u%, one sees that u c a is monotone decreasing. Thus d\u c a < 
and deriving equation (|19p and applying once more the maximum principle gives 
dml < 0. 

It remains to study the behavior of u c a with respect to c. The continuity is 
deduced from the uniqueness of the solution and a priori estimates in the standard 
way. Now let c\ < c 2 and denote by u\ (resp. by U2) the solution of (fT§)) with 
c = c\ (resp. c = 01). Since d\U\ < 0, 

Alii + c 2 diiii + f(ui) - ag(y)ui = (c 2 - ci)d\Ui < 

and U\ > is a super-solution of equation (|19[) with c = C2. By uniqueness of the 
solution, necessarily u 2 < Mi- Once more the strong maximum principle implies 

U 2 < Ml. 

4.2 Convergence to a solution on M. N 

Now that the equation is solved on a domain bounded in the xi-direction, the idea 
is to increase a up to infinity so that the domain tends to M. N . However if c is 
chosen arbitrarily, the function u c a may converge toward the constant or to V 
when a tends to infinity. Hence we adopt a normalization method as in [3] ■ The 
following theorem will define the value of the speed c depending on a to avoid those 
situations. We recall that since a < a>o, A Q the principal eigenvalue of the linearized 
operator about the solution is negative. 

Theorem 4.3. Let us fix e > 0. Let S > be such 
that S < — X a < f'(0). Let 77 > be such that 
Vs £ [O.ry] f(s) > (/' (0) - 6)s. We fix 6 e (0, f). 
Then there exists A(e) > such that for all a > 
A(s), there exists a unique speed c a € (0, 2y/—X a +e] 
with <"(0,0) = 6. 

Proof. By continuity and monotonicity, it suffices to prove: 

i) «°(0,0) >9, 

ii) u°(0, 0) < 6 for c > 2V— Ac + e and a large enough. 

i) Assume c = 0. Let ip a be the positive eigenfunction of the linearized operator 
L associated with the first eigenvalue X a < and with the normalization ||y>a||oo = 
1. Let us introduce v(x\,y) = h(x\)(p a (y) for (xi,y) E [—a, a] x M. N where h(xi) = 
r) 2 ^ 1 . Then < v < 7] on [-a, a] x R N , which yields 

-Av + ag(y)v - f(v) < -Av + ag(y)v - (/'(0) - 5)v 
= (A a + S)v < 0. 

Moreover by construction of V (cf section[3]), v(— a, y) = r]ip a (y) < V(y) if rj is small 
enough. Then v(a, y) = and u < 1. Hence, w is a sub-solution of (|19l) for c = 0. 
Thus u < and therefore, u°(0,0) > u(0,0) = | > 6>. 

ii) Let us construct an explicit super-solution for c > 2y / —X a +e. We recall from 
section T3.ll that Xr < X a for (3 < a and that lim^Q, A^j = A Q . Thus there exists 
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(3 e (0, a) such that 2^—Xp < 2^/—X a + e. As before, let tpp denote the positive 
eigenfunction of the linearized operator L associated with the first eigenvalue A^j < 
with the normalization ||V^i||oo = 1. Choose R such that fro all r > R (a — (3)g(r) + 
Xp > 0, ag(r) > /'(0), and k > is such that kijip > V on Br. The constant k 
only depends on /3 hence on e. 

Lemma 4.4. Then ktpp > V on R N ~ 1 . 

Proof of the lemma. We follow a similar proof to that of lemma 01 
If the lemma does not stand, since ktpp — V tends to at oo, there exists yo <E 
jjiV-i \y -q r such that (hpp - V)(y ) = min R Jv-i (kipp - V) < 0. At this point, 
A(ktpp - V) > but 

-A(ki>p -V) + ag(y)(kiPp - V) - /'(())(*% -V) = 

((a - /%(|t/ |) + Xp)i>p + f'(0)V ~ f(V) > ((a - f3)g(\y \) + Xp)ipp. 
By the choice of R, we get A(kipp — V)(yo) < which yields a contradiction. □ 



Let us now build the super-solution when c > 2y/—X a + e. We set w(x\,y) 
z(xi)kt/)p(y) where z is the solution of 

{z" + cz' — Xpz = on (—a, a), 
z(-a) = l, z(a) = 0. 

Then w verifies 

{-Aw — c<9iu; + ag(y)w = (a — (3)g(y)w + f'(0)w > f(w) 
w(—a, .) = fc?A/3 > tw(o, •) = 0. 

so it is indeed a super-solution of (fP9|) . Moreover 



e p + (a:-a) _ e p-(x-a) 
Z(X) = 5 ^ > 



where p_ < p + < are the roots of p 2 + cp — Xp — 0, i.e. p± = C± ^2 + — (note 
that c 2 + AXp > 0). Hence 

o<<(o,o) < Ko,o) = e e _7 + + 2 I:^L fc#(o) 

e -p + a _|_ g-p-O ^ v y — 

Thus if a is large enough to have e~i a < |, then for any c > 2y/—X a + e, we get 

«S(o,o)<0. □ 

With the bounds on the speed c a it is now possible to pass to the limit as a 
tends to infinity. 

Proposition 5. There exists a sequence (aj)j £ jq such that Oj — > +oo, c aj — > c* G 
[0, 2\J— X a ] and u a " J — > it in C 2 oc (R w ). TTie iimi< n is solution of 

j -An - c*diu + aff(j/)u = / (n) on R w 
\o<u<^, u(O,O) = 0, 9in<0. 

TTien u is necessarily a traveling front solution of (|18[) wif/i c = c* . 
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Proof. First for j G N, fix £j = and aj > A(£j) with — > +00. Since c aj G 

(0, 1\J— \ a + £j] is bounded, u a a / is uniformly bounded in C 2,7 for any 7 G (0, 1). 
Hence up to an extraction of a subsequence, there exist c* G [0, 2\J— \ a \ and u G C 2 oc 
such that c a — > c* and u c a a — > it. Clearly the function u is a solution of ([23]). Owing 
to the normalization u(0, 0) = 8, u is not a constant, moreover by the maximum 
principle < u < V and d x u < 0. Since u is decreasing in x, u± — lim u(x, •) are 

x— >±oo 

solutions of ([10]) and u_(0) > 8 > and < u + (0) < 6. This implies that u_ = V 
and «,+ = 0. Thus u is indeed a traveling front solution of f|18[). □ 



4.3 Existence of traveling front for c > 2\/— A Q . 

In this section we still assume < a < a$ and we will prove the following theorem. 

Theorem 4.5. There exists a traveling front of speed c of equation (|18p if and only 
if c> 2a/~ A q . 

We start with the Proposition 



Proposition 6. For c < 2\/—\ a there exists no traveling front solution of (|18|) . 
Thus c* = 2 v 7 — A a (where c* is the traveling speed constructed in the previous 
section). 

Proof. We argue by contradiction and assume that there exists a traveling front 
u of speed c < 2\/— A Q of (|T8|) . We are going to construct a small positive sub- 
solution with compact support. To this end, we can find 5 G (0,/'(0)) such that 
c 2 + 4(A Q + 28) < and n > such that for all s G [0, 17], f(s) > (/'(0) - <5)s. 

Since the linearized operator L = —A + ag(y) — f'(0) is self adjoint, the principal 
eigenvalue X a is the limit of the Dirichlet principal eigenvalue in Br when R — > 00 
(see [13] for more details): 

f-A^ + a 5 (y)^-f(0)^ = A^, y G Br, 
\if; R >0 onB R , ip R = on dB R . 

Precisely A„ > A Q and \ R jR '^°°> A Q . In the following, ip R denotes the positive 
cigenfunction with ||^ fl ||oo = 1 and let us fix R sufficiently large so that A Q < A^ < 
K + S. 

Let a = a + i-£r, L > 0, be an imaginary root of X 2 + cX — X R — 8 = which is 
possible since c 2 + 4(Af + 8) <c 2 + 4(A Q + 28) < 0. Finally let us fix e > small 
enough such that ee axi < r] and ee tTXl ip R (y) < u(xi,y) for x G [-L, L] and y G Br. 
We set 

fse™i cos(&a?i)^ H (y) if - L < a?i < £, y G 
w(xi,y) = < . . ^ (25) 

otherwise. 



Then w verifies 



-Aw; — c9iw + ag(y)w = (/'(0) — <5)u> < /(w) 



since < w < r\. Moreover w < u and w > on (—L,L) x Br. Thus w is a 
generalized sub-solution with compact support 



17 



Let us now derive a contradiction with the existence of a traveling front u. Trans- 
late u to the left by defining u T [x\, y) = u(x\-\-t, y) for r > 0. Since u(x\, .) X1 ^ r+00 y 
0, there exists r* > such that u T * > w but u T *{x\,y*) = w(x*, y*). Since u T * > 0, 
x* G (-L, L) and y* G Br (an interior point of the support of w). Now since w is 
a sub-solution, the strong maximum principle yields u T * Eicon [— L, L] x but 
this is impossible on the boundary. □ 



We have already proved that for c < 2-^— X a , there exists no traveling front of 
speed c solution of (fl8|) and that for c = c* = 2-^— A a there exists a traveling front 
of speed c. Let us prove that for any c > c* there exists at least a traveling front 
to conclude with Theorem 14.51 The proof goes as usual. We consider the following 
problem 

i-Au-cdiu + ag(y)u = f(u), x = (xi,y) <E (—a, a) x R w_1 ^ 
[u(— a, •) = u*(— a + r, •), u(a, •) = u*(a + r, •) 

where u* is the traveling front of speed c*. The function u*(- +r, •) is a strict super- 
solution of ((26)) (since c > c*) when is a strict sub-solution and < + r, ■). 
Hence as in theorem 14. 1[ it can be proved that there exists a unique solution of 
PB"|) and moreover c^w^ < and 

V(xi ! y) e [-a, a] x M^ 1 V(j/) > u*(-a + r,y) > v r a (x u y) > u*(a + r,y) > 0. 

By uniqueness, w r a depends continuously on r 6 1, sojbJ > and w r a > 

V uniformly on [—a, a] x R^" 1 . Let us denote u a = v r a where r is chosen in order 
that 1^(0,0) = 8 (see previous section for definition of 6). Once again taking any 
sequence a n — > +oo, up to an extraction u an — > u in C^ oc and u is a traveling front 
of speed c solution of (fTS)) . 



5 The case of a KPP non-linearity. Asymptotic 
speed of spreading. 

This section is concerned with the asymptotic behavior of the solutions of the 
parabolic problem 



d t u -Au= f(u) - ag(y)u on R x R N 
u(0, x) = uq(x) on 



(27) 



where / is KPP and uo is an initial condition at least bounded. 
5.1 Extinction for a > a Q 

Let us fix a > ato- We recall that there is no positive asymptotic profile of (fit))) . 

Theorem 5.1. For uq € L°° , there exists a unique solution u(t,x) of (|27p and it 

verifies u(t, x) t ^ +co - > q uniformly for x £ M. N . 

This section is devoted to the proof of this theorem. 
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Let us fix S = max(l, || wo|| oo)- Then the constant functions and 5" are respec- 
tively sub- and super-solutions of (|27|) . Thus there exists u(t,x) a solution of (|27|) 
such that < u < S. By the parabolic maximum principle, this solution is unique. 

Let us define w the solution of (|27|l with the initial condition w(0, x) = S. Since 
the problem and the initial condition do not depend on x\, neither does w thus we 
will write w(t,y). By the maximum principle, < u < w < S and since S is a 
super-solution, dtw < 0. Thus w(t,y) t ^ + °°> W(y) and 

< limsupw < W < S. 

t— >+oo 

Now by parabolic local estimates, W is necessarily solution of 

-A y W = f(W) - ag(y)W 

and thus is a nonnegative asymptotic profile. Since a > ao, W = 0. So u(t, x) 
tends to for t — > +oo uniformly in M. N . 



5.2 Spreading for a < ao 

In this section we assume a < ao- So there exists a critical speed c* of existence 
of traveling front for (JTHJ). We assume that u$ G Cq(M w ), i.e. u is continuous and 
compactly supported, and that uq < V where V is the positive asymptotic profile 
solution of (|TU)) . We will prove the spreading of the solution of (f^7|) but we first 
need the following theorem. 

Theorem 5.2. The unique solution of 

j-Az-cd 1 z + ag(y)z = f(z) (x 1: y)e™ N 
}0< z(xi,y) <V(y) (x u y) e 



(28) 



with c < c* is z(x,y) = V(y). 

Proof. Let us consider the generalized sub-solution with compact support w(x\, y) 
defined in (|2"5]l . This is possible since c < c*. Up to a decrease of e > 0, we can 
assume that w < z on M. N . Now by applying the sliding method to w T where 
w T (xi 1 y) = w(xi + T,y) and z, one can prove that w T < z for all r£l. We can 
thus define 

My e K"" 1 z(y) = inf z{x lt y) > 

and state that z ^ 0. Now z is a super-solution of (fT0|) since z = inih^s. z(- + h, •) 
and an infimum of solutions is a super-solution. 

Finally as in section [3J we can build a positive sub-solution of (fTUf smaller than 
z and thus by monotone iteration we have a solution of ([10[) between these sub- 
and super-solution. By uniqueness of the positive solution, we obtain V < z. And 
due to condition in (|28|) . we have z = V. □ 



Let us now turn to the precise study of the spreading of the solution of (|27j) 

Theorem 5.3. For u £ Cq(M. n ) with u < V, there exists a unique solution u of 
((27)) and 

for any c > c* lim sup u(t, x) = 0, (29) 

*^+°° |xi|>ct 

for any c with < c < c* lim sup \u(t,x) — V(y)\ = 0. (30) 

t ^+ 00 |x 1 |<ct 
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Proof. Fix c > c* . Let U denote a traveling front of speed c* . Since U{x\ , •) — > V for 
2:1 — > —00 locally uniformly, there exists L £ R such that f/(xi — L,y) > iio(xi,y) 
for all (xi, y) € l w . Now considering v(i, x) = U{x\ — L — c*t 7 y) and applying the 
comparison principle, we have u(t, x) < i>(t, x) for alH > and x € R w . Thus since 
U is decreasing in xi 

sup u(t, x) < sup U(x\ — L — c*t,y) = sup U{{c — c*)t — L,y). 

\xi\>ct \xi\>ct 



Since c > c* and U(xi,y) Jl ^ + °°> uniformly in y € We see that 

supj. > ct x) — ► as t — > +00. Since u(t, —xi,y) satisfies the same equation 
(|2"7|) . this shows that sup a , 1< _ ct u(t, x) — >■ as well as t — > +00. Thus (j2"i?j) is proved. 
Assume now c < c* . Let us first prove the following weaken version of (|30p : 



Lemma 5.4. For any c € R with |c| < c*, 

V(xi,y)eR JV lim Ki,xi-ct,2/)-F(2/)| =0. (31) 

£— ► +00 

Proof of lemma \5^4\ Let us assume that c > 0, the proof being similar for c < 0. 
Let ?;(£, x±,y) = u(t, x\ — ct, y). Then v satisfies the equation 

d t v — Av — cd\v + ag(y)v = f(v) (32) 

with the initial datum v(0,xi,y) = uo(x,y) > and ^ 0. Hence by the parabolic 
maximum principle, for all (xi,y) £ R w v(l,xi,y) > 0. Now since c < c* , in 
(|25|) . we constructed w(xi,y) > a stationary non-zero sub-solution of (|32|) with 
compact support and w could be chosen arbitrary small. Hence we can assume 
w 5: f(l, •, •). So if w is the solution of 

{dtw — Aw — cdiw + ag(y)w = f(w) t > 0, (xi,y) £ R w 
w(0,xi,y) = w(xi,y) (x l7 y)£R N 

then by comparison principle, Vt > 1 V(xi, y) £ M. N v(t, x±, y) > w(t— 1, xi, y). Now 
since w is a sub-solution, w is increasing with respect to t and < w{t,X\,y) < 
V(y). Therefore, by standard elliptic estimates, w(t,x\,y) t ~ s,+00 > z(x,y) and z is 
a solution of (|28p . By theorem 15.21 we have z = V and this complete the proof 
of the lemma since by the comparison principle w(t — l,x\,y) < v(t,xi,y) < V(y) 
thus 

lim v(t,xi,y) = V{y) 
which yields ((311) • □ 

Let us now prove (|30|) , that is the uniform convergence to V in the expanding slab 
{xi < ct}. We will only prove it for < x\ < ct. Indeed using as before u(t, — Xi, y), 
the general result follows from the convergence in the set {0 < xi < ct}. 

Let c with < c < c* be fixed and let e > be given (arbitrarily small) . For R > 
sufficiently large, we know that the principal eigenvalue A^f of the problem (f24|) 
above is such that < 0. Denote by tp R > the corresponding eigenfunction of 
(|24]l. Under these conditions we know that there exists a unique solution V (y) > 
of the profile equation in Br with Dirichlct condition: 

f -AV R + ag(y)V R - f(V R ) in B R 

\V R : = on V R > in 1 J 
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(Compare e.g. [J]). Moreover, it is straightforward to show that V R is increasing 
with R and that limfl_j. +00 V R (y) = V(y). 

Let us choose R > sufficiently large so that for all y Br V(y) < £ and for 
all y E Br < V(y) — V R {y) < e. The proof of the uniform convergence to V for 
c < c* will rest on the following Proposition. 

Proposition 7. Let c be such that < c < c* . Then, with R chosen as above, there 
exists a solution v c (x\, y) defined for x\ £ R~ , y € Br of equation 

— Av — cd\v + ag(y)v = f(v) x\ < 0, y € Br (34) 

satisfying the following properties: 

' v c > and d\v c < in R~ x Br, 
v c (0,y) = for y € Br, 
v c (xi,y) =0fory<E dB R , xi < 0, 
v c (-oo,y) = V R (y) for y g B R . 

Postponing the proof of this proposition, let us complete the proof of Theorem 
Extending v c by for x\ > turns v c into a (generalized) sub-solution of 
equation (f3"4"]) in the cylinder R x Br (see [TJ). Therefore u c (a;i — c(t — t ),y) is a 
sub-solution of the equation (f27|) in this cylinder for all to > and all c £ (0, c*). 

By Lemma [5.41 (applied here in the case c = 0), we can fix to > sufficiently 
large such that for t > t we have 

u(t,0,y)>V(y)-- for all ?/ g ~B~r 

where S = min (V - V R ) > 0. Therefore, 

Br 

u(t, 0, y) > V R (y) for all i > t and all y G Sr. 

We fix c g (c, c*) and we consider x\,y) = Wg(^i —c(t — to),y). In the region 
_D = (0, +oo) x Br, u is a solution and v a sub-solution of equation ([27} and for 
any time t > to 

u(t, xi,y) > v(t, x\,y) for (xi,y) e dD. 
Moreover, u{to,x\,y) > v{to,x\,y) = in D. The comparison principle then yields 

u{t,x 1 ,y)>v{t,xi,y) in D. 

Therefore 

limsup sup (^(y) — u(i, a;)) < limsup sup (V(y) — Vz(x\ — c(t — t ), y)) 

< limsup sup [Vc((c— c)t + cto,y)) 

t-t+oo y£B R 

< sup (V(y) - V R {y)) < e. 
yeB R 

Outside of Br we already know that < u < V < e for any t > 0, x\ € R and 
\y\ > R. Therefore 

limsup sup (V(y) — u(t, x)) <e 

t->+oo 0<X!<ct 
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ee axi cosi^X! + if - 2L < x x < 0, y G B R , 

otherwise. 



Since this is true for all e > (and for — ct < x\ < 0), we have thereby established 

It now remains to prove Proposition [7] which we carry now. As in (|25[) , we 
construct a sub-solution of the equation (|34j) with compact support, namely: 

In comparison with (I25p . there is a translation in xi such that the support of w 
now lies in E_ x Br. 

For any b < 0, let z b be the solution of 

-Az b - cdiz b + ag(y)z b = f(z b ) in (6, 0) x B R , 
z b {b,y) = V R {y), z fc (0,y) = for y G TT^ 
Z b (x 1 ,y) = for x\ G (b, 0), \y\ = R. 

Since y fl is a super-solution and a sub-solution, there exists a solution of this 
problem. By the sliding method of [3], we know that this solution is unique and 
satisfies d\z b < in (6, 0) x Br. 

Next, for b < —L, wince w is a sub-solution, we also know that 

Vb<-L V(xi,y) E (b, 0)xB ^(xi, y) > w{x u y). 

This allows us to pass to the limit when b — > — oo. Clearly 0t(a;i, y) — ^> w c (^i, y)- 
By the lower bound, v c (x\, y) > w(xi, y) which shows that v c (x\, y) > in Ml x Br. 
Since 9i« c < and v c ^ 0, we also know that 9iu c < in ML x _B^. Now since 
Iinx^-^-oo v c (xi, y) must be a positive solution of (|33j) . Hence by uniqueness we 
get v c (—oo,y) = V' R (j/). This completes the proof of Proposition [7] and therefore of 
Theorem [Ql □ 



6 The case of a bistable non-linearity. Asymptotic 
profiles. 

In this section we consider again equation (TTJ) but in the bistable framework. That 
is, we assume that / is a C 1 function that satisfies the following assumptions for 
some 6 G (0, 1): 

/(0) = f(0) = /(l) = 0, f(s) < for a G (0, 0) and f(s) > in (0, 1), (35) 

/'(0) > 0, f'(l) > 0. (36) 

We also assume that 

f f(s)ds>0. (37) 
Jo 

We are concerned here with the existence of traveling front solutions of ([1]), that is, 
(c, u) solution of (fT8|) . First we require some preliminary results on the equation 
(fT0|) in the bistable case. 
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6.1 Existence of asymptotic profiles in the bistable case 

Consider equation 

Au + f(u)-ag(y)u = 0, y€R N ~\ 
u > 0, u bounded, 

under the same assumption ([7]) and ([8} as above for the function g. 

The existence of solutions depends on a and is obtained in the following theorem. 

Theorem 6.1. Let f and g satisfy the above assumptions. There exists a threshold 
value a* £ (0,oo), such that: 

i) For any a G (a*,+oo), (|38p does not have any positive (non-zero) solution. 

ii) For any a G (0, a*] } (|38[) admits a maximal positive solution V(y). 

iii) For any a G (0, a*), (|38|) admits a second positive solution W(y) with < 

W(y) < V(y). 

The rest of this section is devoted to the proof of this Theorem. 

This Theorem follows from the observation that for a > any positive solution 
u(y) of Q38p satisfies u(y) — > as \y\ — > oo. This is obtained from Corollary[TJ 

Next, by the maximum principle, any solution of (|38[) satisfies < u < 1 (we 
think of f(s) as having been extended by outside [0, 1]). 

Now u = 1 is a super-solution of problem ([38| . Any solution of (|38| for a 
is a sub-solution of (|3"B"|) for any parameter /3 < a. Therefore, if there exists a 
positive bounded solution of (|38p for a, there also exists a positive solution for any 
< P < a. 

Next, we claim that for small enough a > 0, (|38j) admits a positive solution. 
Indeed, consider the functional defined on H: 

J(w) = J a (w) = [ (l\Vw\ 2 + ^g(y)w 2 - F(w)) dy 

where F(z) = J* f(s)ds. Recall that / is extended by outside [0, 1], thus F is 
bounded. Since g(r) ~ > oo as r —> oo, it is straightforward to show that there exists 
a minimizer v of J(w): J(v) — min{J(u>), w G %}. Furthermore, we know that 
v > and v is a solution of (|38| (see Theorem 13.31 for details). 

Let us show that for a > small enough J(v) < 0. To this end, let be defined 

by 

f l if \y\ < R 

(n(y) = <R + l-\y\ ifR<\y\<R + l 
[o ii\y\>R+l 

Then ( R G H and 

MCr) = I - HCr) < -f(i)\b r \ + C\B R+1 \ B R \ 

where \A\ denotes the volume of A and C is a constant. Since — F(l) < by ([37)) . 
we see that by choosing R large enough, Jo(Cr) < 0. Then for such an R fixed, 



23 



we see that JaiCn) < provided a > is small enough. This guarantees that 
J a (v) < 0. 

It follows that v ^ 0. By the maximum principle, we then have < v < 1. This 
shows that for small a > 0, (|38| admits a positive solution. 



Next, we show that if a is large enough (|38[) does not admit any positive solution. 
This can be seen by multiplying the equation by u and integrating to yield: 



|Vup + a / ff (2/K = / f( u ) u <mju z (39) 
f(s) 

where m = sup > 0. We conclude with the following lemma: 

s>0 S 

Lemma 6.2. Under the assumption ([8]) g(r) r ~ i ' 00 ) oo, for any e > 0, £/iere exists 
a constant K(e) > smc/i i/iai /or all u one /ias: 

/ u 2 <e[ \S7u\ 2 + K{e)[ g{y)u 2 . 

Jr n ~ 1 Ji"- 1 it"- 1 

Indeed choosing in the lemma e = we get from (|39]l 

i | |Vn| 2 + L- m if(^-)) Jg(y)u 2 <0 

This shows that for a > mif (^-), the only solution of (f38|) is u = 0. 

Proof of Lemma lKM Let <5 = (5(e) > be chosen such that the principal eigenvalue 
of —A in Hq(P>2s) is larger than |. Let % be a smooth cutoff function such that 
x{r) = 1 if < r < S, xi r ) = if r > 25 and < x < 1. Consider ui = \ u an d 
U2 = (1 — x) u so that u = u\ +U2- Using (a + b) 2 < 2 (a 2 + b 2 ) , since Ui € Hq(E>2s) 
by Poincare's inequality, we have 

m 2 =j u x 2 < e -\ |v Ul | 2 <ff/ |V«|V+/ « a |Vxl a Y 

JB 2i 4 iB 2i 1 \Jb 2I JB 2S \B s ) 

So that 

,2 



J R N-1 I J B J\, 



u 

\v\>s 
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where ki(e) > |Vx| 

Next / u 2 2 < u 2 . Therefore 



j i u 2 < 2 (j Ul 2 + J u 2 2 ^j < e J \Vu\ 2 + K{e) J g{y)u 2 
where K(e) = 2 efc y^ +1 . The lemma is thus proved. □ 



The next step is to prove that the set of a > such that (|38j) has a solution is 
a closed set. Let atj — >• a* be a sequence such that (|38|) admits a solution itj such 
that < < 1 for all j. Note that by the maximum principle, 9 < maxiij < 1. 
The sequence (iij) is bounded by 1 and by standard elliptic estimates is locally 



24 



compact. Therefore, one can extract a subsequence Uj such that Uj — > u* uniformly 
on compact sets in the C 2 -norm. Therefore u* is a solution of (|38[) for the value 
a = a*. We know that u* > 0, but since max itj > 0, we see that max it* > 9. 
Indeed by section [21 Uj{y) — > as \y\ — > oo uniformly with respect to j. Therefore 
u* > and (j38|) also has a positive solution for a* . This shows that the set of a 
such that (p?5|) has a positive solution is an interval (0, a*] with < a* < oo. 
Considering the evolution equation 

d t z-Az + ag(y)z = f(z), t>0,yeR N ~\ 
z(0,y) = l, 

we see that t h> z(t,y) > is decreasing and therefore has a limit. This limit 
is necessarily the maximal positive solution V — V a for the a for which ([551) has 
a positive solution, that is a S (0,a*], or is in the opposite case, that is when 
a > a . 

The existence of a second solution when < a < a* is inspired from a work 
of P. Rabinowitz [M] . In a slightly different formulation, the existence of pairs of 
solutions is established in [33] by a topological degree argument for bistable type 
nonlinearities and another type of parameter dependance. The use of the topological 
degree involves compact operators and the results of [34] are set in the framework 
of bounded domains. A similar construction can be carried here owing to the 
condition ([8]) g(r) — > +oo as r — > +oo. Indeed, under this condition, the injection 
W L 2 (R N ^ X ) is compact. This allows one to construct a compact operator and 
to carry the argument of [34j to the present framework. 

Since we will not use the second solution, we will leave out the details of the 
proof of the existence of a second solution. 

6.2 Stable asymptotic profiles 

As we have seen, a solution of (|38[) is obtained by the minimization of J = J a 
defined above. The proof of the existence of the previous solution for a > small 
yields the following result. 

Proposition 8. There exists < a* < a* such that for all a € (0,a*) there exists 
a minimum v a > of J a and such that 

Ja(va) = rnin J a < 0. 

In the following, we require the notion of stable solution. 

Definition 6.3. Let v be a solution of ([38jl . Eigenvalues of the linearized problem 
about v are defined as the eigenvalues A of 

-A(p + ag(y)(p- f(v)(p = \<p in M. N ~ l . 

The principle eigenvalue is uniquely determined by the existence of a corresponding 
eigenfunction (p with <p > 0. We say that v is (weakly) stable if the principal 
eigenvalue A = Xi[v] of the linearized problem satisfies Xi[v] > 0. 

It is well known that the maximal solution V(y) given by Theorem 16.11 when 
< a < a* is weakly stable. Likewise, the minimum solution of the energy of the 
Proposition [8j above, when < a < a*, is a weakly stable solution. 
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In the following we consider the case < a < a* and we make the following 
assumption. 

There exists a unique positive stable solution of (fTU|) . (40) 

This condition implies that the minimizer solution v a : J a {v a ) = min{ J a (v), v G TL\ 
coincides with the maximum solution V. 

We leave it as an open problem to give sufficient conditions for the uniqueness 
of the stable solution. Uniqueness results have been given for analogous problems 
but with a = 0, which would rather correspond to the minimal solution in our 
framework |32j . Likewise it would be interesting to give sufficient conditions that 
ensure that a* = a*. Condition (j40|) has several implications that wc can state. 

Proposition 9. For a £ (0, a*], under condition f|40j) . there does not exist a pair 
of distinct ordered functions (vi,v 2 ) with < V\ < Vi < V, v\ is a sub-solution 
and V2 is a non-maximal solution. That is, if < V\ < v 2 < V are respectively 
sub-solution and solution of (|38p . then v\ = v 2 - 

Proof. The proof follows the observation in However, it requires new elements 
in view of the unbounded domain. If v\ < V2, let ip2 be a principal eigenfunction of 
the linearized problem corresponding to Ai[u2]. Since and V are the only stable 
solutions, Ai[«2] < 0. We claim that for e > sufficiently small, v = v% — eip2 is a 
super-solution of (|55|) . Indeed 

-Av + ag(y)v - f(v) = f{v 2 ) -/(«)- /' '{v 2 )eip 2 - Ai [v 2 ]eip 2 

( f(v 2 ) - f(v 2 - eip 2 ) ,,, . . r A 

= / (V2) - Ai [V 2 \ £ip 2 - 

V eL Pl ) 

The right hand side is positive if e > is sufficiently small. 

Next, given R > 0, we can choose e > small enough so that v\ < V2 — stp 2 in 
Br. We choose R so that v\{y) < 5 for all \y\ > R and / is decreasing on [0,5]. 
We claim that then v\ < V2 — eip2 in R^ 1 \ Br. Argue by contradiction. In this 
were not the case, then, since v\, V2 and if2 converge to at infinity, there exists y, 
\y\ > R such that 

min{w 2 -e(p 2 -vi} = v 2 {y) - Sip 2 (y) - Vi(y) < 

R"-l 

This implies that < v(y) < v±(y) < 5. Denote L the operator L = —A + ag(y). 
Since < L(v - u x ) - (f (v) - /(ui)) and f(v(yj) - f{vi(y)) > 0, at the point y 
we get L(v — vi)(y) > 0. Therefore, we have reached a contradiction. This shows 
that vi < V2 — sif2- Now we have a super-solution v above a sub-solution v\. This 
implies that there exists a stable solution v such that v% < v < V2 — £<f2 < V. This 
however is in contradiction with condition (|40j) . □ 

From this property, we derive the following useful consequence. 

Proposition 10. Let a £ (0, a*) and let W be the maximal solution of equation 
(|38[) with the value a* of the parameter. Then, any other solution v of (|38[) with 
parameter a that is not the maximal solution cannot be above W . 

This immediately follows from the previous proposition as W is a sub-solution 
of the equation for the value a < a* and W < V. 
A consequence of this proposition is 
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Proposition 11. For a £ (0,a*] and under condition (|40[) . the maximal solution 
V is isolated in L°° topology. Therefore, there exists 9\ > 9 such that if v is a 
solution of (|38p with v(0) > 9\ then v = V . 

As we have done before, we can prove that if v is a solution such that v > W in 
Br, then v > W < V in R w_1 . Then any solution v ^ V is such that there exists 
y € -Br such that v(y) < W(y) and therefore \\v — > min-g— V — W = 5 > 0. 

Now, if there exist a sequence v n of solutions of ([55)1 such that t> n (0) — ► V(0) 
then by elliptic estimates v n — > W a positive solution of ([55)1 and W(0) = V^(0) so 
W = V by the maximum principle which contradicts the fact that V is isolated. 



7 Traveling fronts for a bistable non-linearity 



In this section we assume that / if of bistable type and satisfies (|35 |) - (|37|) . In 
addition, we assume that < a < a* and that condition (|40| is fulfilled. Therefore, 
there exists a unique non-zero stable solution V(y) = V a (y) of the profile equation 
([55]). Therefore V > 0, J a (V) = min{J(», w G %}, Ai[V] > and V" is isolated 
in the L°° topology. Furthermore 1^ is the maximal solution. Any other non-zero 
solution w satisfies < w < V in R^" 1 and Ai[w] < where Ai[u>] is the principal 
eigenvalue of the linearized problem defined in definition 16.31 

In this section, we prove the existence of a traveling front solution of (TTJ) repre- 
senting an invasion of by the state V at positive speed. Such a solution is given 
as a pair (c, u) of 



— Ait — cd\u + ag{y)u = f(u) in 
u(-oo, y) = V(y), u(+oo, y) = 



(41) 



with c < and u : R N -> (0, 1). 

We follow the construction of a solution given above. Namely, let a > 1 and in 
the slab E a = (—a, a) x R" -1 , consider the problem 

j - Am - ediu + ag(y)u = f(u) in E , 
|u(-a,y) = y(y), u(+a,j/) = 0. 

We recall that for any c g R, for a fixed, there exists a unique solution u = u c 
of (|42[) . Furthermore, < u < V and d\u < in S a .The mapping c t— > w c is 
decreasing. 

Up to here, the procedure is the same as before. From this point on however, 
we need to modify the above argument since we used the fact that / was positive. 
Our first task is to prove the following 



Proposition 12. There exists a unique (c a ,u a ) such that u a is a solution of (|42[) 
for speed c a and u a satisfies the normalization condition 

max u a (0,y) = 9. (43) 

Let us first prove the existence of c a . The uniqueness is clear. 

The parameter c a is bounded independently of a > 1 . Moreover, 

liminf c a > 0. 

a— > +oo 
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Proof. The bound from above is obtained simply by comparison with the one di- 
mensional problem. Indeed, consider the ODE problem for z = z(x\): 



| - z " - 1*' = /0) in (~a, °) f44) 
\z(-a) = l, z(+a) = 0, z{0) = 9 

It is known that there exists a unique value j a for which (|4"4")l has a (unique) solution 
z. Furthermore, lim a _j. +00 7° = 7* where 7* is the unique speed of traveling fronts 
for the ID equation 

{— z" — 7*2' = f(z) in (—a, a) 
z(-oo) = 1, z(+oo) = 

Comparing (|44p with (|4"2"j) , we see that for each = 7°, the solution z of ([4"4"]) is a 
super-solution of (gU), thus z > u 7 ° and for all y G R N_1 , it 7 °(0,j/) < z(0) = 0. 
Since c 1— >■ tt c is decreasing, we see that 

max u c (0, y) < 8 for all c > y a . 

Assume now that max. y(Efl jv-i u c (0,y) < 8 for all c E R. Passing to the limit for 
c — > —00, u c converges toward a positive solution v of with maxu < 6*. By the 
maximum principle, it is impossible thus there exists a unique c a £ (—00,7°) such 
that gSl) is fulfilled. 

Since 7° — > 7* < 00 as a — > +00 and a t— > 7° is a continuous function, this shows 
that sup c a < 00 . 

a>l 

Since a 1— > c a is continuous, in order to complete the proof of the Proposition, 
it suffices to show that liminf a _>. 00 c a > 0. For this, we argue by contradiction and 
assume that for a sequence aj — > +00 there holds c aj < 0. For the sake of simplicity, 
we write a instead of aj. Since c 1— > u c is decreasing, from this we infer that along 
this subsequence, the solution v = v a of 

{—Av + ag(y)v = f(v) in S a 
v{-a,y) =V(y), v(+a,y)=0 

satisfies max ygR jv-i v(0,y) < 8. 

Due to Proposition [Til there exist 8\ > 8 such that if an asymptotic profile v 
solution of (f38|) verifies v(0) > 8\ then v = V. 

There is a point b — bj, —a < b < such that v a (b, 0) — Q\. We now translate 
the solution to center it on X\ = b. That is, we let v a (xi, y) — v a (xi + b, y) defined 
for x\ G (—a — b, a — b) and y E The interval (—a — b,a — b) either converges 

(along a subsequence) to (—00, +00) or to some (— d, +00) with < d < 00. In both 
cases, by standard elliptic estimates, one can strike out a subsequence of v a , denoted 
again v a , such that v a converges locally to some function w where w satisfies: 

i-Aw + ag(y)w = f(w) in (~d, +00) x R w_1 
[diiv^O, w(0,0)=9 1 . 

In case the interval is converging to (— d, +00), in addition we know that w(—d, y) = 
V(y). If the interval converges to M, then lim^-^-oo w(x\, y) exists and is some 
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function W(y) which is then a solution of the profile equation (|38[) . But since 
w(0,0) = Oi, we know that W(0) > 9\. By the definition of 6\, this implies that 
W = V . Therefore, denoting d = oo in case (—a — b, a — b) — > R, in both cases, we 
get 

VyeM^" 1 w{-d,y)=V(y) 

where now < d < +00. We also know that w(+oo, y) = ip(y) exists with < ip> < 
V. 

Multiply by d\w and integrate over (— d, +00) x to get 

f l(diw) 2 + Jty)- J(V) =0 

J{x 1= -d] * 

where diw = if d = 00. In all cases, we get 

j(v) > 

Since V minimize J a , we obtain V = ip and w(a;i, y) = V(y) for all x\ G (— d, +00) 
but this contradicts the renormalization u>(0, 0) = 0\. 

Wc have thus reached a contradiction. This shows that for large a, c a > 0, which 
completes the proof of the Proposition. □ 

Let us now turn to the proof of the existence of traveling front solutions of 
(|4T|) . Since c a and u a are bounded, by standard elliptic estimates, we can strike 
out a sequence a = a,j — ¥ 00 (we continue to denote subsequences by a) such that 
c a — > c > and u a — > u. We know that (c, u) satisfies the equation 

—Am + cd\u + ag(y)u = f(u) in ~R N 

with d\u < and max fi »-i u(0, •) = 9. It remains to identify the limits as x\ — > 
±00. These ]im x - L -^± oc u(x\, y) = u±(y) exist and are solutions of the asymptotic 
profile equation ((55)) . Now since < u + (y) = \im xi ^ +00 u(xi,y) < 9 and all 
positive solutions w of satisfy m&xw > 9, we have u + = 0. 

We claim that it_(y) = lim ;ci _ > ._ 00 u(x\, y) coincides with V(y). Clearly, < 
ii_ < V. Argue by contradiction that u_ ^ V, implying U— < V. By assumption, 
U— is an unstable solution of (|38|) in the sense that Ai[u_] < 0. Let us construct a 
super-solution of the stationary equation, that is a w with 

—Aw + ag{y)w > /(u>) 

such that to is a compact perturbation of it_ and as close as we wish to u~ . 
Consider the linearized equation about it_ : 

-Aip - f'(u-(y))tp + ag(y)4> = \x[u-]ip 

with Ai[u_] < 0. We know that Ai[u_] is the limit of the Dirichlet principal 
eigenvalue in a ball when the radius goes to infinity (This follows from the Rayleigh 
quotient minimization). Therefore, R > can be chosen sufficiently large so that 
the principal eigenvalue [i and associated cigenfunction ip of 

-Alp + ag(y)ip - f'(u-)ip = flip in B fl 
ip = on dB R , V > in Br 
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satisfy fx < 0. 

Consider the function £(xi,y) = cos(u>xi)ijj(y) defined for x\ G (—L,L) with 
L = £j and \y\ < R. We note D = (—L,L) x Br. This function is positive and 
satisfies: 

-AC + ag(y)C- f(u^C = ( f , + ^K in D 
C = on 3D 

Choose L large enough so that /i + uj 2 < 0. Then let w(x\,y) = U-{y) — e((xi,y) 
with e > and (xi,y) G D. This function satisfies 

-Aw + ag{y)w - f(w) = (-(ft + co 2 ) + /("-) ~ /("- ~ ^0 fl[u ^ ^ 

Since fi + uj 2 < 0, we can choose e sufficiently small so that 

—Aw + ag(y)w — f(w) > in D and w > 0. 

Furthermore, because e( = on dD and e£ > in D, that is w < u_ in D, if we 
extend w by choosing w(xi,y) = U-(y) for all (xi,y) ^ £*, we have constructed a 
(generalized) super-solution of the problem (see e.g. [TTj). 
Let us now derive a contradiction. We consider two cases. 

• Cases (i): Suppose c > 0. Then U(t,X\,y) = m(xi — ct,y) is a solution of the 
evolution equation 

d t U - AU + ag(y)U = f(U) t G R, (xi.y) e R*. 

Now [/ *~ 5 ' °°> locally uniformly in (xi , y) . Furthermore, for all times U(t,xi,y) < 
u-(y). Since w is a compact perturbation of u_ for a time to sufficiently negative, 
we get 

V(xi,y)€R N U(t ,x 1 ,y)<w(x 1 ,t). 

Now when U(t,xi,y) t ^ + °°y u-(y) locally uniformly and we get a contradiction 
since U(t, x±,y) < w(xi,y) < U-(y) for all (x±,y) G D. 

• Case (ii) : The case that remains to be studied is c = (since we already have 
c > 0). Then u(x±, y) is a stationary solution of the same equation that w is a super- 
solution of. Since u(— oo, y) = u_(y) and u(+oo, y) = 0, and since w = u_ outside a 
compact set, after a translation, we can assume that = u(x\+h, y) < iy(xi, j/)(for 
large enough h). Define 

h* = m£{h G R,u(x! + h,y) < w{xx,y) in R N }. 

Clearly, h* > — oo (for to < it_ at some points). Then w{x\, y) > w(xi + /i*, j/) = u^* 
and min(ui — u^*) — is necessarily achieved at a point of D. Since w(xi,y) = 
U-(y) > u(x\ + h, y) for all h if (xi, y) $ D, we see that the maximum is achieved 
at an interior point of D. Writing w — u^* > as a super-solution of a linear elliptic 
equation in D, we derive a contradiction with the strong maximum principle. 
Therefore in all cases, the solution u satisfies the limiting condition: 

w(-oo, y) = V(y), u(+oo, y) = 0. 

Therefore (c, u) is a solution of the traveling front equation (|4T]) . 
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8 The model of cortical spreading depression 



We consider here more general versions of the model ([5]) described in the Introduc- 
tion. The problems studied in this paper have the following general form 

d t u- Au = h(y,u) x=(x u y)£R N . (46) 

In the modeling context N = 2 and 3 are the cases of interest. As indicated in 
the Introduction, this equation also describes cortical spreading depressions (CSD). 
There the wave propagates in a medium composed of two different components, the 
gray and white matters of the brain, with a narrow transition area separating them. 
Thus we consider in this section functions h(y, u) of the following type: 

h(y,u)=f(u)fov\y\<L 1 (47) 

h(y^ u ) < —mu for \y\ > L 2 (48) 

I y I — ^^oo _i_ 
h(y, u) + mu — > uniformly for u G M (49) 

where < L\ < L 2 < oo and K > m > are given parameters and / is of bistable 
form. That is we assume that / verifies conditions (|3"5j) - ([3"T)) of section [5] Note that 
in particular, we assume 

/ f(s)ds>0. 
Jo 

We also assume that y >-> h(y, s) is continuous and that s >-> h(y, s) is Lipschitz 
continuous for all s <G [0, 1] (and \y\ ^ L\ in case L\ = L 2 ). Lastly we assume that 

Vs g [0, 1] Vy g R w_1 h(y, s) < max{/(s), -ms}. 



8.1 The asymptotic profile equation 

We start as usual with the profile equation 

(-AV = h(y,V) yeR N -\ 
\V>0, Abounded. 

We recall that Ai[V] is the principal eigenvalue of the linearized equation about V. 
This eigenvalue can be defined as 

i / ?6-ff 1 (R JV ~ 1 ) J (fi 2 

Associated with ([50| is the energy functional: 



J(w) 



[ (l\Vw\ 2 - H(y,w)) dy 



where H(y,z) = h(y,s)ds. Note that owing to condition (|48|). J(w) is well 
defined for all w £ H 1 (R N ^ 1 ). 

Theorem 8.1. There exist critical radii < < L* < oo with the following 
properties: 
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i) For L/2 < L*, there is no solution other than to the asymptotic profile equation 
Mb- 



ii) For L\ > L* (independently of L2), there exists a maximum solution V of (|50[) 

and this solution is stable in the sense that Xi[V] > 0. 

iii) For all L\ > L* , the minimum of the energy functional is achieved at some 

non- zero function Vj(y), i.e. 

JiVj) = min J(w) < 0. 

weH 1 (R N - 1 ) 

Proof, i) Since the equation implies that — Au + mu < for all \y\ > L2 in R w_1 , 
and u > is bounded, by Theorem 12. II we know that u and \Vu\ have exponential 
decay as \y\ — > +00. Then we get 



a 2 . 



min(l, m)||u||iji(RJv-i) < / |Vu| 2 +mw 2 < / (f(u)-\-mu)u<K 

Jr"- 1 Jb L2 

We know, by Sobolev embedding and Holder inequality, that 

u 2 < n(L 2 )\\u\\ hi^n-i) 

where 77(^2) — > as L2 —> 0. Therefore for L2 small enough, these inequalities yield 
u = 0. 

ii) Next, since 1 is a super-solution of the equation in M^ -1 , there exists a 
maximum solution of equation ([50]) that we denote V. By what we have just seen, 
V = for L2 sufficiently small. Let us now show that V > for L\ sufficiently 
large. 

Let us consider the energy restricted to the ball of radius R < L\ 



Jr{w) 



where F(z) = J Q f(s)ds. We know (see the proof of Theorem 16. lj) that for R 
sufficiently large there exists a minimum u>r of 

Jr(wr) = inf Jr{w) < 0. 

Then wr > in Br and wr is solution of —Avjr = /(wr) in Br and wr = on 
8Br. Extending wr by outside the ball Br, we get a global (generalized) sub- 
solution. The solution V is such that V > wr (since V is the maximum solution). 
This implies that V ^ and therefore V > in R N - 1 for Li > R. 
iii) Now for L\ > R, clearly 

inf J(w) < Jr(wr) < 0. 
Let us now show that the infimum is achieved. 



32 



Let (w n ) be a minimizing sequence: J(w n ) — > inf J < for n — > +00. Note that 
J is bounded from below. Writing 



J{wn) > 0/ |Vw„| 2 + mw„ 2 + / i|Vw„| 2 - H (y, w n ) 

z Jr n - 1 \b L2 Jb l . 2 1 

If \S7w n \ 2 + mw„ 2 + [ i|Vw„| 2 - C + ew n 2 

1 JR«-1\-Br., JBr,„ Z 



> - 

~ 2 



'\Bl 

= -C\B L2 \ + e||w„|| ff i( R jv-i) 

we can strike out a subsequence still denoted (w n ) such that w n — > w weakly in 
H 1 (R N ~ 1 ). Now using (gHJ), for every e > 0, there exists R = R(e) > such that 
\F(y, s) + y s2 \ < £s2 f° r a ^ \v\ — R ana ^ au s £ ^ + (there is no loss in generality 
in assuming w n > as w n + is also a minimizing sequence). Therefore 

J(wn) = x I |Vu; n | 2 + ^ / w„ 2 + r(e) - / H(w n , y) 

Z JK N -! Z J\y\>R J\v\<R 

where \r(e)\ < Ce for some constant C > 0. 

By compact injection of i? 1 (R Ar_1 ) L 2 (Br), we can assumme that w n — > w 
strongly in L 2 (Bl 2 ). Then by standard arguments relying on Lebesgue's dominated 
convergence Theorem, we see that 



H(y,w n ) / H(y,w). 

Br J Br 

Owing to Sobolev embedding and H older inequality, it is straightforward to check 
that the quantity 

\ [ |V W „| 2 + ^ f w n 2 
z Jr™- 1 z J\y\>R 

defines the square of a norm equivalent to the usual _ff 1 (R JV_1 ) norm. Hence using 
the lower semi-continuity of the norm, we get 

lim J(w„) = inf J > - / |Vw| 2 + ^ / w 2 +r(e)-[ H(y,w). 

n ^°° 2 J R N-1 2 J\y\> R J Br 



Now using again (|49[) . we get: 



JH= / |V W | 2 - / H(y,w)< inf J + 2|r(e)|. 

Since e > is arbitrarily small we get JYuO = inf J. □ 

Remark 2. By using the method of [2], one can show that for L\ large enough 
there exists a second solution in R^ -1 . 



We will now make use of the condition that the stable solution of (|50|) is unique. 
In the paper of Chapuisat and Joly |18) . it is argued by phase plane method, that 
for the case N — 1 = 2, L\ = L 2 and f(y, s) = —ms for |y| > L 2 that indeed this is 
the case. We note that it is an interesting open problem to derive such uniqueness 
results in more general situations or to complete the heuristic part of the argument 
of E51. 
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8.2 Traveling fronts for the CSD model 



In this section, we prove Theorcm ll.4l The proof is similar as in Scction[7] There we 
used that h(y, it) = /(it) — ag{y)u with g — > +00. But actually, the same properties 
that were entailed one can derived for h(y,u) < — mu for large \y\. We start by 
constructing a solution of 

{Au a -c a d 1 u a = h(y,u a ) in (-a,a) x l^ 1 
\u a {-a, y) = V(y), u a {+a, y) = 0, 

that verifies 

sup u a (Q,y) = 9 (52) 

for a > 1 and where 9 is the unstable of /(it), that is f(0) = and < 9 < 1. Wc 
recall that c a is uniquely determined by the renormalization condition (1521) . 

< I S NJ- /(«) 

Let /(it) = max{/(u), — m(it)}. Note that / itself is / \ 

bistable: _\ „/ \, 



1(0) = f(9) = /(l), f(s) < in (0, 9), f(s) > in (9, 1). 
We denote by z c a the solution of 

{-z"-cz> = f(z) 
\z(-a) = l, z(+o) = 0. 



(53) 



The function z^ a is a supcrsolution of (|5"Tj) thus u a < z^ a . In view of (f52|) this 
implies that z° a (0) > and this implies that c a < 7a where 7 a is the unique value 
of c such that the solution of (|53|) verifies z^" (0) = 0. This as before yields the 
upper bound for c a . 

The lower bound is achieved in the same manner as in the section [7] and the 
convergence for a — > +00 also. 
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